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IONOSPHERIC  MODELING 

E  and  F  Region  (100-500  km) 

A  computer  code  has  been  developed  that  will  calculate  the  time- 
varying,  vertical  ionization  distribution  between  100  and  500  kilometers 
altitude  for  the  polar  ionosphere.  The  code  produces  ionization  from 
both  the  solar  ultraviolet  flux  and  by  auroral  electron  precipitation 
down  the  geomagnetic  field  lines.  The  distribution  of  ionization  into 
the  various  ion  species  is  traced  through  ion-chemical  reactions,  which 
are  specified  by  the  user.  At  present,  66  reactions  can  be  accommodated. 
The  ion-chemical  reactions  form  a  set  of  production  and  loss  terms,  which 
can  be  folded  into  the  species  continuity  equations  to  calculate  the 
time  variable  species  density.  A  typical  continuity  equation  can  be 
written 

3N4 

5T*  pi  -  L,  -  »<*<«, V  (1) 

where  P^,  L.. ,  V.. ,  and  N..  are  total  production  rate,  loss  rate,  velocity 
and  number  density  for  the  1th  species.  An  equation  such  as  (1)  must  be 
written  for  each  species  whose  density  Is  varying  with  time;  a  complete 
solution  for  n  species  requires  the  simultaneous  solution  of  n  continuity 
equations.  The  divergence  term  In  equation  (1)  represents  the  rate  of 
change  In  the  local  density  due  to  motion  of  particles.  This  term  may 
take  on  a  variety  of  forms,  depending  upon  the  complexity  of  completeness 
one  wishes  to  specify  the  species  densities  and  their  velocities.  If 
motion  Is  Ignor'd,  Dlv(N^V^)  may  be  considered  zero,  and  the  resulting 
equation  Is  a  first-order  ordinary  differential  equation.  If  molecular 
diffusion  Is  desired,  then  the  divergence  term  takes  the  font 
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3N .  32N> 

Dtv(N.V.)  -  a(2)Ni  +  b(2)^  +  dzj-jpj-  (2) 

where  a,  b,  and  c  are  coefficients  that  depend  upon  the  kinetic  pro¬ 
perties  of  the  gas;  l.e. ,  temperature,  velocity,  mobility,  collision 
rate,  etc.,  and  gravity.  Substituting  (2)  into  (1)  produces  a  second- 
order  partial  differential  equation  which  is  much  more  difficult  to  solve 
than  an  ordinary  differential  equation.  If  three-dimensional  winds  and 
pressure  gradients  are  included  in  the  model,  the  divergence  operator 
produces  partial  derivatives  in  three  coordinate  directions,  and  the 
resultant  partial  differential  equation  in  3  space  variables  and  1 
time  variable  becomes  extremely  difficult  and  time  consuming  to  solve 
numerically.  Our  Initial  calculations  neglect  all  motion.  Thus,  these 
calculations  will  be  of  limited  use  except  to  evaluate  production  and 
loss  terms  and  to  provide  Ionization  profiles  which  can  be  compared  with 
profiles  obtained  with  vertical  diffusion.  In  order  to  Include  vertical 
diffusion,  the  subroutine  that  solves  the  system  of  continuity  equations 
will  have  to  be  replaced.  A  copy  of  such  a  computer  code  is  now 
in  house  and  is  presently  being  evaluated. 

A  flow  chart  of  the  computer  program  Is  shown  in  figure  1.  The  user 
has  the  option  of  equal  altitude  Increments  or  random  altitudes.  Densities 
for  0,  0,,  and  H,  are  generated  from  the  atmospheric  model  of  Jacchia  [1971]. 
Initial  densities  for  the  other  species  in  the  system  may  be  input  or 
are  assumed  to  be  10"4  times  the  H,  density.  Data  on  the  atom  composition 
and  total  concentr  lion  and  on  charged  species  are  input  so  that  checks 
on  the  charge  and  i/ass  balance  can  be  made. 

Ionization  i;  p^odjeed  by  both  solar  flux  and  auroral  particles. 

These  processes  ft  further  subdivided  into  daytime  and  nighttime 
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solar  fluxes  and  into  drizzle  particles  and  arc  particles. 

The  calculation  of  ionization  rates  by  solar  ultraviolet  radiation 
is  based  on  the  derivation  originally  presented  by  Chapman  [1931]  and 
extended  to  a  multispecies  gas  by  Swider  [1964].  To  calculate  the  ionization 
rate  at  a  given  altitude,  it  is  necessary  to  integrate  the  absorption 
of  the  solar  flux  from  the  top  of  the  atmosphere  down  the  ray  path  to 
the  height  of  interest.  The  net  solar  flux  that  remains  is  then  available 
to  ionize  whatever  neutral  atmospheric  species  that  might  exist  at  the 
given  altitude.  The  difficulty  with  this  calculation  is  that  the  solar 
spectrum  is  composed  of  both  emission  bands  and  lines  and  each  wavelength 
interval  must  be  handled  separately.  In  addition,  the  solar  absorption 
and  ionization  cross-sections  of  the  atmospheric  gases  show  considerable 
structure  and  it  is  necessary  to  carefully  correlate  cross-section  data 
in  order  to  match  the  cross-sections  to  the  solar  lines  and  bands.  Above 
100  km  the  photoionization  and  photoabsorption  cross-sections  of  Stolarski 
and  Johnson  [1972]  have  been  used,  since. these  cross-sections  correspond 
with  the  solar  flux  lines  and  bands  measured  by  Hinteregger  [1970].  These 
data  cover  the  solar  spectrum  between  30-1027  A,  which  is  the  important 
part  of  the  spectrum  for  the  E  and  F  layers. 

A  nighttime  photoionization  source  has  also  been  Incorporated  into 
the  computer  code.  This  source  Is  based  on  solar  hydrogen  and  helium 
emission  lines,  which  are  resonantly  scattered  to  the  nighttime  Ionosphere 
by  the  hydrogen  and  hcllua  gases  of  the  earth's  upper  atmosphere.  Ogawa 
and  Toheatsu  [1966]  first  suggested  the  brightness  of  these  sources 
(HLya  1026  A--10  raylelghs,  Ulya  1216  A--4  kllorayleighs ,  He!  584  A--10 
raylelghs,  Mel  I  304  A— 1  raylelghs)  and  both  Chen  and  Harris  [1971]  and 
fujitaka  et  al.  [1971]  shoved  tl  it  this  scattered  radiation  would  be 


sufficient  to  account  for  the  preservation  of  the  nighttime  E  layer. 

The  initial  nighttime  model  is  simplified  by  assuming  a  constant  emission 
brightness  throughout  the  night,  but  this  will  eventually  be  replaced  with 
a  calculation  of  the  multiple-photon  scattering  in  spherical  geometry. 

Calculation  of  the  auroral  ionization  rates  is  based  on  an  empirical 
relation  developed  by  Lazarev  [1967]  from  laboratory  data.  It  has  been 
found  experimentally  that  when  a  monoenergetic  beam  of  electrons  is 
directed  to  a  gas,  such  as  N2,  02,  air,  etc.,  that  the  rate  of  ionization 
and  depth  of  penetration  into  the  gas  by  the  beam  depend  on  the  density 
of  the  gas,  the  initial  energy  of  the  electrons  in  the  beam,  and  the 
flux  of  the  electrons,  but  are  insensitive  to  the  composition  of  the  gas. 
It  is  also  found  that  on  the  average,  each  electron  in  the  beam  loses 
about  34  eV  per  collision  as  it  slows  down.  An  expression  was  developed 
by  Lazarev,  based  on  these  observations  and  represents  an  empirical 
fit  to  laboratory  data,  v/hich  calculates  the  ionization  rate  and  depth 
of  penetration  for  an  electron  of  a  given  energy.  To  rake  Lazarev's 
equation  applicable  to  the  Ionosphere,  we  need  to  specify  the  energy 
spectrum  (flux  as  a  function  of  energy)  of  the  Incoming  electrons  at 
the  "top"  of  the  atmosphere.  Rocket  and  satellite  observations  have 
shown  that  the  energy  spectrin  of  auroral  electrons  Is  highly  variable; 
therefore,  the  computer  program  has  been  written  to  accept  any  energy 
spectrum  as  Input  data,  or  to  generate  Its  own,  according  to  the  type  of 
aurora  desired  (ICC  class). 

The  calculation  proceeds  as  follows:  The  electron  flux  as  a  function 
of  energy  Is  cither  read  In  or  generated  Internally  as  an  exponential 
function  of  energy,  for  the  case  of  the  Internal  particle  spcctnaw,  a 
para«cter  mst  be  specified  to  designate  the  IK  type  aurora.  These 


electron  fluxes  are  the  fluxes  at  the  "top"  of  the  atmosphere  (500  km). 

The  mass  density  of  the  atmosphere  is  then  integrated  downward  from 
500  km  and  at  each  height  the  penetration  depth  and  ionization  rate  are 
calculated  from  Lazarev's  equation.  Summing  over  the  entire  energy 
spectrum  yields  the  total  ionization  rate,  which  is  then  proportioned 
between  the  species  in  proportion  to  their  relative  abundances. 

The  low-energy  electron  fluxes  that  produce  the  "drizzle"  ionization 
are  based  on  the  electron  spectrum  observed  by  Burch  [1968]. 

The  computer  code  has  been  written  so  that  the  photoionization  due 
to  the  solar  flux  is  time  dependent  through  the  solar  zenith  angle. 

The  "drizzle"  Ionization  is  presently  assumed  to  be  constant  with  both 
ti  e  and  latitude.  The  auroral  arc  ionization  is  turned  on  and  off  at 
times  specified  by  the  user.  This  allows  simulation  of  an  arc  of  any 
duration  at  any  time  du-ing  the  day  or  night. 

Thus,  ionization  profiles  between  100-500  kra  can  be  calculated  at 
a  given  latitude  throughout  the  day.  By  choosing  different  latitude 
values,  a  complete  ionization  nap  can  be  generated  for  the  polar  ionosphere. 


0  Region  (30-100  Lwl 

In  our  attack  or.  the  0  region  problem,  we  are  following  a  philosophy 
in  which  thc  sophistication  of  the  node!  grows  in  tic*.  As  a  conscience, 
our  current  work  Is  involved  with  the  completion  of  «.ny  of  the  subroutines 
that  we  intend  t0  use.  Three  production  sources  art  planned  for  incorporation 
into  the  0  region  model.  These  Include  modifications  of  the  solar  photos 
absorption  and  the  auroral  tlearen  ca1ulatl*m  as  they  are  no*  being 
developed  for  the  r  region  problt*  in  addition,  «e  er*  edanfnf  projection 
rate  calculations  Are  to  solar  rep***  „4  alp*,*  stifles  as 
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by  Adams  and  Masley  [1965].  At  a  later  time,  it  is  planned  that 
brehmstralung  effects  will  also  be  included;  for  some  types  of  events, 
this  may  be  an  important  low-altitude  ionization  source.  At  this  time 
we  will  also  extend  the  electron  ionization  calculation  to  include 
relativistic  effects. 

Electron  density  calculations  will  involve  the  computation  of  electron 
loss  rates.  In  order  to  be  completely  general,  this  term  must  depend 
upon  time,  altitude, and  electron  density.  The  plan  is  to  begin  the  program 
with  a  semi -empirical  set  of  these  coefficients  calculated  by  using 
experimental  profiles  and  an  extension  of  the  two-ion  model  as  described 
by  Adams  and  Meg ill  [1967]  as  a  means  of  interpolating  between  the 
various  sets  of  data  that  are  available.  The  extension  will  essentially 
create  a  four-ion  model  since  two  types  of  positive  ions  will  be  considered 
in  order  to  take  into  account  the  fast  recombination  rates  of  the  "hydrated 
positive  ions.  Subsequent  versions  of  this  calculation  will  include 
production  of  positive  ion  profiles  for  comparison  with  mass  spectrometric 
data  where  desirable.  It  is  planned  to  have  a  working  program  by  the  end 
of  the  summer  with  future  work  going  toward  the  extensive  ion  chemistry 
calculations  and  sophistication  of  the  production  rate  programs. 

At  present  It  is  not  planned  to  attempt  a  complete  Ion  and  neutral 
ch«*1siry  calculation  in  a  single  progra*.  The  ainor  species  chemistry  In 
the  neutral  cheaistry  models  Is  poorly  known  and,  as  a  consequence,  the 
itiAlure  of  these  nodcls  would  perhaps  confuse  the  issue  rare  than  clarify 
(l.  t*e  plan,  instead,  to  wake  ettenslvc  Hies  of  the  diurnal  and  perhaps 
seasonal  variations  of  the  neutral  species  and  use  these  as  a  basis  for 
0*r  ion  chemistry  calculations. 

Hf.-re  7  it*  computer  program  outline. 


Figure  2.  Computer  Program  Outline  for  the  0  Region 
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